3.484 \(\int (a+a \sin (e+f x)) \sqrt{c+d \sin (e+f x)} \, dx\)

Optimal. Leaf size=179 \[ -\frac{2 a \left (c^2-d^2\right ) \sqrt{\frac{c+d \sin (e+f x)}{c+d}} F\left (\frac{1}{2} \left (e+f x-\frac{\pi }{2}\right )|\frac{2 d}{c+d}\right )}{3 d f \sqrt{c+d \sin (e+f x)}}-\frac{2 a \cos (e+f x) \sqrt{c+d \sin (e+f x)}}{3 f}+\frac{2 a (c+3 d) \sqrt{c+d \sin (e+f x)} E\left (\frac{1}{2} \left (e+f x-\frac{\pi }{2}\right )|\frac{2 d}{c+d}\right )}{3 d f \sqrt{\frac{c+d \sin (e+f x)}{c+d}}} \]

[Out]

(-2*a*Cos[e + f*x]*Sqrt[c + d*Sin[e + f*x]])/(3*f) + (2*a*(c + 3*d)*EllipticE[(e - Pi/2 + f*x)/2, (2*d)/(c + d
)]*Sqrt[c + d*Sin[e + f*x]])/(3*d*f*Sqrt[(c + d*Sin[e + f*x])/(c + d)]) - (2*a*(c^2 - d^2)*EllipticF[(e - Pi/2
+ f*x)/2, (2*d)/(c + d)]*Sqrt[(c + d*Sin[e + f*x])/(c + d)])/(3*d*f*Sqrt[c + d*Sin[e + f*x]])

________________________________________________________________________________________

Rubi [A]  time = 0.197774, antiderivative size = 179, normalized size of antiderivative = 1., number of steps used = 6, number of rules used = 6, integrand size = 25, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.24, Rules used = {2753, 2752, 2663, 2661, 2655, 2653} \[ -\frac{2 a \left (c^2-d^2\right ) \sqrt{\frac{c+d \sin (e+f x)}{c+d}} F\left (\frac{1}{2} \left (e+f x-\frac{\pi }{2}\right )|\frac{2 d}{c+d}\right )}{3 d f \sqrt{c+d \sin (e+f x)}}-\frac{2 a \cos (e+f x) \sqrt{c+d \sin (e+f x)}}{3 f}+\frac{2 a (c+3 d) \sqrt{c+d \sin (e+f x)} E\left (\frac{1}{2} \left (e+f x-\frac{\pi }{2}\right )|\frac{2 d}{c+d}\right )}{3 d f \sqrt{\frac{c+d \sin (e+f x)}{c+d}}} \]

Antiderivative was successfully verified.

[In]

Int[(a + a*Sin[e + f*x])*Sqrt[c + d*Sin[e + f*x]],x]

[Out]

(-2*a*Cos[e + f*x]*Sqrt[c + d*Sin[e + f*x]])/(3*f) + (2*a*(c + 3*d)*EllipticE[(e - Pi/2 + f*x)/2, (2*d)/(c + d
)]*Sqrt[c + d*Sin[e + f*x]])/(3*d*f*Sqrt[(c + d*Sin[e + f*x])/(c + d)]) - (2*a*(c^2 - d^2)*EllipticF[(e - Pi/2
+ f*x)/2, (2*d)/(c + d)]*Sqrt[(c + d*Sin[e + f*x])/(c + d)])/(3*d*f*Sqrt[c + d*Sin[e + f*x]])

Rule 2753

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)]), x_Symbol] :> -Simp[(d
*Cos[e + f*x]*(a + b*Sin[e + f*x])^m)/(f*(m + 1)), x] + Dist[1/(m + 1), Int[(a + b*Sin[e + f*x])^(m - 1)*Simp[
b*d*m + a*c*(m + 1) + (a*d*m + b*c*(m + 1))*Sin[e + f*x], x], x], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*
c - a*d, 0] && NeQ[a^2 - b^2, 0] && GtQ[m, 0] && IntegerQ[2*m]

Rule 2752

Int[((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])/Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]], x_Symbol] :> Dist[(b*c
- a*d)/b, Int[1/Sqrt[a + b*Sin[e + f*x]], x], x] + Dist[d/b, Int[Sqrt[a + b*Sin[e + f*x]], x], x] /; FreeQ[{a
, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && NeQ[a^2 - b^2, 0]

Rule 2663

Int[1/Sqrt[(a_) + (b_.)*sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Dist[Sqrt[(a + b*Sin[c + d*x])/(a + b)]/Sqrt[a
+ b*Sin[c + d*x]], Int[1/Sqrt[a/(a + b) + (b*Sin[c + d*x])/(a + b)], x], x] /; FreeQ[{a, b, c, d}, x] && NeQ[a
^2 - b^2, 0] &&  !GtQ[a + b, 0]

Rule 2661

Int[1/Sqrt[(a_) + (b_.)*sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[(2*EllipticF[(1*(c - Pi/2 + d*x))/2, (2*b)
/(a + b)])/(d*Sqrt[a + b]), x] /; FreeQ[{a, b, c, d}, x] && NeQ[a^2 - b^2, 0] && GtQ[a + b, 0]

Rule 2655

Int[Sqrt[(a_) + (b_.)*sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Dist[Sqrt[a + b*Sin[c + d*x]]/Sqrt[(a + b*Sin[c +
d*x])/(a + b)], Int[Sqrt[a/(a + b) + (b*Sin[c + d*x])/(a + b)], x], x] /; FreeQ[{a, b, c, d}, x] && NeQ[a^2 -
b^2, 0] &&  !GtQ[a + b, 0]

Rule 2653

Int[Sqrt[(a_) + (b_.)*sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[(2*Sqrt[a + b]*EllipticE[(1*(c - Pi/2 + d*x)
)/2, (2*b)/(a + b)])/d, x] /; FreeQ[{a, b, c, d}, x] && NeQ[a^2 - b^2, 0] && GtQ[a + b, 0]

Rubi steps

\begin{align*} \int (a+a \sin (e+f x)) \sqrt{c+d \sin (e+f x)} \, dx &=-\frac{2 a \cos (e+f x) \sqrt{c+d \sin (e+f x)}}{3 f}+\frac{2}{3} \int \frac{\frac{1}{2} a (3 c+d)+\frac{1}{2} a (c+3 d) \sin (e+f x)}{\sqrt{c+d \sin (e+f x)}} \, dx\\ &=-\frac{2 a \cos (e+f x) \sqrt{c+d \sin (e+f x)}}{3 f}+\frac{(a (c+3 d)) \int \sqrt{c+d \sin (e+f x)} \, dx}{3 d}-\frac{\left (a \left (c^2-d^2\right )\right ) \int \frac{1}{\sqrt{c+d \sin (e+f x)}} \, dx}{3 d}\\ &=-\frac{2 a \cos (e+f x) \sqrt{c+d \sin (e+f x)}}{3 f}+\frac{\left (a (c+3 d) \sqrt{c+d \sin (e+f x)}\right ) \int \sqrt{\frac{c}{c+d}+\frac{d \sin (e+f x)}{c+d}} \, dx}{3 d \sqrt{\frac{c+d \sin (e+f x)}{c+d}}}-\frac{\left (a \left (c^2-d^2\right ) \sqrt{\frac{c+d \sin (e+f x)}{c+d}}\right ) \int \frac{1}{\sqrt{\frac{c}{c+d}+\frac{d \sin (e+f x)}{c+d}}} \, dx}{3 d \sqrt{c+d \sin (e+f x)}}\\ &=-\frac{2 a \cos (e+f x) \sqrt{c+d \sin (e+f x)}}{3 f}+\frac{2 a (c+3 d) E\left (\frac{1}{2} \left (e-\frac{\pi }{2}+f x\right )|\frac{2 d}{c+d}\right ) \sqrt{c+d \sin (e+f x)}}{3 d f \sqrt{\frac{c+d \sin (e+f x)}{c+d}}}-\frac{2 a \left (c^2-d^2\right ) F\left (\frac{1}{2} \left (e-\frac{\pi }{2}+f x\right )|\frac{2 d}{c+d}\right ) \sqrt{\frac{c+d \sin (e+f x)}{c+d}}}{3 d f \sqrt{c+d \sin (e+f x)}}\\ \end{align*}

Mathematica [C]  time = 6.28371, size = 1736, normalized size = 9.7 \[ \text{result too large to display} \]

Warning: Unable to verify antiderivative.

[In]

Integrate[(a + a*Sin[e + f*x])*Sqrt[c + d*Sin[e + f*x]],x]

[Out]

a*((c*Sec[e]*(1 + Sin[e + f*x])*(-((AppellF1[-1/2, -1/2, -1/2, 1/2, -((Csc[e]*(c + d*Cos[f*x - ArcTan[Cot[e]]]
*Sqrt[1 + Cot[e]^2]*Sin[e]))/(d*Sqrt[1 + Cot[e]^2]*(1 - (c*Csc[e])/(d*Sqrt[1 + Cot[e]^2])))), -((Csc[e]*(c + d
*Cos[f*x - ArcTan[Cot[e]]]*Sqrt[1 + Cot[e]^2]*Sin[e]))/(d*Sqrt[1 + Cot[e]^2]*(-1 - (c*Csc[e])/(d*Sqrt[1 + Cot[
e]^2]))))]*Cot[e]*Sin[f*x - ArcTan[Cot[e]]])/(Sqrt[1 + Cot[e]^2]*Sqrt[(d*Sqrt[1 + Cot[e]^2] + d*Cos[f*x - ArcT
an[Cot[e]]]*Sqrt[1 + Cot[e]^2])/(d*Sqrt[1 + Cot[e]^2] - c*Csc[e])]*Sqrt[(d*Sqrt[1 + Cot[e]^2] - d*Cos[f*x - Ar
cTan[Cot[e]]]*Sqrt[1 + Cot[e]^2])/(d*Sqrt[1 + Cot[e]^2] + c*Csc[e])]*Sqrt[c + d*Cos[f*x - ArcTan[Cot[e]]]*Sqrt
[1 + Cot[e]^2]*Sin[e]])) - ((2*d*Sin[e]*(c + d*Cos[f*x - ArcTan[Cot[e]]]*Sqrt[1 + Cot[e]^2]*Sin[e]))/(d^2*Cos[
e]^2 + d^2*Sin[e]^2) - (Cot[e]*Sin[f*x - ArcTan[Cot[e]]])/Sqrt[1 + Cot[e]^2])/Sqrt[c + d*Cos[f*x - ArcTan[Cot[
e]]]*Sqrt[1 + Cot[e]^2]*Sin[e]]))/(3*f*(Cos[e/2 + (f*x)/2] + Sin[e/2 + (f*x)/2])^2) + (d*Sec[e]*(1 + Sin[e + f
*x])*(-((AppellF1[-1/2, -1/2, -1/2, 1/2, -((Csc[e]*(c + d*Cos[f*x - ArcTan[Cot[e]]]*Sqrt[1 + Cot[e]^2]*Sin[e])
)/(d*Sqrt[1 + Cot[e]^2]*(1 - (c*Csc[e])/(d*Sqrt[1 + Cot[e]^2])))), -((Csc[e]*(c + d*Cos[f*x - ArcTan[Cot[e]]]*
Sqrt[1 + Cot[e]^2]*Sin[e]))/(d*Sqrt[1 + Cot[e]^2]*(-1 - (c*Csc[e])/(d*Sqrt[1 + Cot[e]^2]))))]*Cot[e]*Sin[f*x -
ArcTan[Cot[e]]])/(Sqrt[1 + Cot[e]^2]*Sqrt[(d*Sqrt[1 + Cot[e]^2] + d*Cos[f*x - ArcTan[Cot[e]]]*Sqrt[1 + Cot[e]
^2])/(d*Sqrt[1 + Cot[e]^2] - c*Csc[e])]*Sqrt[(d*Sqrt[1 + Cot[e]^2] - d*Cos[f*x - ArcTan[Cot[e]]]*Sqrt[1 + Cot[
e]^2])/(d*Sqrt[1 + Cot[e]^2] + c*Csc[e])]*Sqrt[c + d*Cos[f*x - ArcTan[Cot[e]]]*Sqrt[1 + Cot[e]^2]*Sin[e]])) -
((2*d*Sin[e]*(c + d*Cos[f*x - ArcTan[Cot[e]]]*Sqrt[1 + Cot[e]^2]*Sin[e]))/(d^2*Cos[e]^2 + d^2*Sin[e]^2) - (Cot
[e]*Sin[f*x - ArcTan[Cot[e]]])/Sqrt[1 + Cot[e]^2])/Sqrt[c + d*Cos[f*x - ArcTan[Cot[e]]]*Sqrt[1 + Cot[e]^2]*Sin
[e]]))/(f*(Cos[e/2 + (f*x)/2] + Sin[e/2 + (f*x)/2])^2) + ((1 + Sin[e + f*x])*Sqrt[c + d*Sin[e + f*x]]*((-2*Cos
[e]*Cos[f*x])/(3*f) + (2*Sin[e]*Sin[f*x])/(3*f) + (2*(c + 3*d)*Tan[e])/(3*d*f)))/(Cos[e/2 + (f*x)/2] + Sin[e/2
+ (f*x)/2])^2 + (2*AppellF1[1/2, 1/2, 1/2, 3/2, -((Sec[e]*(c + d*Cos[e]*Sin[f*x + ArcTan[Tan[e]]]*Sqrt[1 + Ta
n[e]^2]))/(d*Sqrt[1 + Tan[e]^2]*(1 - (c*Sec[e])/(d*Sqrt[1 + Tan[e]^2])))), -((Sec[e]*(c + d*Cos[e]*Sin[f*x + A
rcTan[Tan[e]]]*Sqrt[1 + Tan[e]^2]))/(d*Sqrt[1 + Tan[e]^2]*(-1 - (c*Sec[e])/(d*Sqrt[1 + Tan[e]^2]))))]*Sec[e]*S
ec[f*x + ArcTan[Tan[e]]]*(1 + Sin[e + f*x])*Sqrt[(d*Sqrt[1 + Tan[e]^2] - d*Sin[f*x + ArcTan[Tan[e]]]*Sqrt[1 +
Tan[e]^2])/(c*Sec[e] + d*Sqrt[1 + Tan[e]^2])]*Sqrt[(d*Sqrt[1 + Tan[e]^2] + d*Sin[f*x + ArcTan[Tan[e]]]*Sqrt[1
+ Tan[e]^2])/(-(c*Sec[e]) + d*Sqrt[1 + Tan[e]^2])]*Sqrt[c + d*Cos[e]*Sin[f*x + ArcTan[Tan[e]]]*Sqrt[1 + Tan[e]
^2]])/(3*f*(Cos[e/2 + (f*x)/2] + Sin[e/2 + (f*x)/2])^2*Sqrt[1 + Tan[e]^2]) + (2*c*AppellF1[1/2, 1/2, 1/2, 3/2,
-((Sec[e]*(c + d*Cos[e]*Sin[f*x + ArcTan[Tan[e]]]*Sqrt[1 + Tan[e]^2]))/(d*Sqrt[1 + Tan[e]^2]*(1 - (c*Sec[e])/
(d*Sqrt[1 + Tan[e]^2])))), -((Sec[e]*(c + d*Cos[e]*Sin[f*x + ArcTan[Tan[e]]]*Sqrt[1 + Tan[e]^2]))/(d*Sqrt[1 +
Tan[e]^2]*(-1 - (c*Sec[e])/(d*Sqrt[1 + Tan[e]^2]))))]*Sec[e]*Sec[f*x + ArcTan[Tan[e]]]*(1 + Sin[e + f*x])*Sqrt
[(d*Sqrt[1 + Tan[e]^2] - d*Sin[f*x + ArcTan[Tan[e]]]*Sqrt[1 + Tan[e]^2])/(c*Sec[e] + d*Sqrt[1 + Tan[e]^2])]*Sq
rt[(d*Sqrt[1 + Tan[e]^2] + d*Sin[f*x + ArcTan[Tan[e]]]*Sqrt[1 + Tan[e]^2])/(-(c*Sec[e]) + d*Sqrt[1 + Tan[e]^2]
)]*Sqrt[c + d*Cos[e]*Sin[f*x + ArcTan[Tan[e]]]*Sqrt[1 + Tan[e]^2]])/(d*f*(Cos[e/2 + (f*x)/2] + Sin[e/2 + (f*x)
/2])^2*Sqrt[1 + Tan[e]^2]))

________________________________________________________________________________________

Maple [B]  time = 0.967, size = 657, normalized size = 3.7 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((a+a*sin(f*x+e))*(c+d*sin(f*x+e))^(1/2),x)

[Out]

2/3*a*(4*((c+d*sin(f*x+e))/(c-d))^(1/2)*(-(-1+sin(f*x+e))*d/(c+d))^(1/2)*(-d*(1+sin(f*x+e))/(c-d))^(1/2)*Ellip
ticF(((c+d*sin(f*x+e))/(c-d))^(1/2),((c-d)/(c+d))^(1/2))*c^2*d-4*((c+d*sin(f*x+e))/(c-d))^(1/2)*(-(-1+sin(f*x+
e))*d/(c+d))^(1/2)*(-d*(1+sin(f*x+e))/(c-d))^(1/2)*EllipticF(((c+d*sin(f*x+e))/(c-d))^(1/2),((c-d)/(c+d))^(1/2
))*d^3-((c+d*sin(f*x+e))/(c-d))^(1/2)*(-(-1+sin(f*x+e))*d/(c+d))^(1/2)*(-d*(1+sin(f*x+e))/(c-d))^(1/2)*Ellipti
cE(((c+d*sin(f*x+e))/(c-d))^(1/2),((c-d)/(c+d))^(1/2))*c^3-3*((c+d*sin(f*x+e))/(c-d))^(1/2)*(-(-1+sin(f*x+e))*
d/(c+d))^(1/2)*(-d*(1+sin(f*x+e))/(c-d))^(1/2)*EllipticE(((c+d*sin(f*x+e))/(c-d))^(1/2),((c-d)/(c+d))^(1/2))*c
^2*d+((c+d*sin(f*x+e))/(c-d))^(1/2)*(-(-1+sin(f*x+e))*d/(c+d))^(1/2)*(-d*(1+sin(f*x+e))/(c-d))^(1/2)*EllipticE
(((c+d*sin(f*x+e))/(c-d))^(1/2),((c-d)/(c+d))^(1/2))*c*d^2+3*((c+d*sin(f*x+e))/(c-d))^(1/2)*(-(-1+sin(f*x+e))*
d/(c+d))^(1/2)*(-d*(1+sin(f*x+e))/(c-d))^(1/2)*EllipticE(((c+d*sin(f*x+e))/(c-d))^(1/2),((c-d)/(c+d))^(1/2))*d
^3+d^3*sin(f*x+e)^3+c*d^2*sin(f*x+e)^2-d^3*sin(f*x+e)-c*d^2)/d^2/cos(f*x+e)/(c+d*sin(f*x+e))^(1/2)/f

________________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int{\left (a \sin \left (f x + e\right ) + a\right )} \sqrt{d \sin \left (f x + e\right ) + c}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+a*sin(f*x+e))*(c+d*sin(f*x+e))^(1/2),x, algorithm="maxima")

[Out]

integrate((a*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c), x)

________________________________________________________________________________________

Fricas [F]  time = 0., size = 0, normalized size = 0. \begin{align*}{\rm integral}\left ({\left (a \sin \left (f x + e\right ) + a\right )} \sqrt{d \sin \left (f x + e\right ) + c}, x\right ) \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+a*sin(f*x+e))*(c+d*sin(f*x+e))^(1/2),x, algorithm="fricas")

[Out]

integral((a*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c), x)

________________________________________________________________________________________

Sympy [F]  time = 0., size = 0, normalized size = 0. \begin{align*} a \left (\int \sqrt{c + d \sin{\left (e + f x \right )}} \sin{\left (e + f x \right )}\, dx + \int \sqrt{c + d \sin{\left (e + f x \right )}}\, dx\right ) \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+a*sin(f*x+e))*(c+d*sin(f*x+e))**(1/2),x)

[Out]

a*(Integral(sqrt(c + d*sin(e + f*x))*sin(e + f*x), x) + Integral(sqrt(c + d*sin(e + f*x)), x))

________________________________________________________________________________________

Giac [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int{\left (a \sin \left (f x + e\right ) + a\right )} \sqrt{d \sin \left (f x + e\right ) + c}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+a*sin(f*x+e))*(c+d*sin(f*x+e))^(1/2),x, algorithm="giac")

[Out]

integrate((a*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c), x)